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A Computational Tool to Interpret the Bulk Composition of Solid Exoplanets 
based on Mass and Radius Measurements 
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ABSTRACT 



The prospects for finding transiting exoplanets in the range of a few to 20 M® is 
growing rapidly with both ground-based and spaced-based efforts. We describe a pub- 
lically availalble computer code to compute and quantify the compositional ambiguities 
for differentiated solid exoplanets with a measured mass and radius, including the mass 
and radius uncertainties. 



Introduction 



Over 250 extrasolar planets are known to orbit nearby main sequence stars. Among these in- 
clude over a dozen exoplanets with minimum masses below 22 and several with minimum masses 
less than 10 Af®. Of key interest are transiting planets with measured masses and radii, which can 
be used to constrain the planet's interior bulk composition. The relationship bet ween mass and 



radius for solid exoplanets has hence received much attention in the last few ye ars ([Valencia et al 



20061 : iFortnev et al 



20071 : ISeaeer et al 



2007 



Selsis et al 



2007 



Sotin et al 



ity builds on much earlier work ( Zapolskv &: Salpeterl 1969 ; Stevenson 1982 ) , with improvements on 



20071 ). The recent activ- 



the equations of state and treatment of different mantle and core compositions to varying degrees 
of complexity. 

Unlike for the solar system planets, we have no access to the gravitational moments of exo- 
planets. Hence the density distribution in the interior is unknown and this leads to an ambiguity, 
or degeneracy, in the interior composition for an exoplanet of a fixed mass and radius. One way to 
capture the degeneracies of exoplanet interior composition is using ternary diagrams (introduced 
to exoplanet interiors by I Valencia et al.1 (|2007l )). 



We adopt the idea of using ternary diagrams to quantify the compositional uncertainty in exo- 
planets. The planet mass and radius are the observed quantities and therefore we foc us solely on on 
ternary diagrams for a planet of fixed mass and fixed radius (c.f. IValencia et al.ll2007l ). We compute 
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ternary diagrams for solid exoplanets ranging in mass from 0.5 to 20 M®. We furthermore explain 
the behavior of the mass-radius curves in two and three dimensions. We also present a description 
of our publically available computer code to compute fixed mass-radius ternary diagrams, including 
the observational uncertainties. 



2. Computer Model 

2.1. Background and Equations 

We begin by assuming the major components of a solid exoplanet are limited to an iron core, 
a silicate mantle, and a water ice outer layer. In other words, we assume the interior of the planet 
is differentiated with the denser materials interior to the less dense materials. We further assume 
each layer to be homogeneous in its composition. 

We can then model the interior of a solid exoplanet by using: 

(1) the equation for mass of a spherical shell 

^ = (1) 

where m(r) is the mass included in radius r, p(r) is the density at radius r; 

(2) the equation for hydrostatic equilibrium 

dP(r) —Gm(r)p(r) 
dr r 2 

where P(r) is the pressure at radius r; and 

(3) the equation of state (EOS) that relates P and p. 

The EOS is different for each different material. We used Fe (e) for the planet core, MgSiOs 
perovskite fo r the s ilicate mantle, and water ice VII, VIII, and X for the water-ice outer layer. See 



Seager et al.l (120071 ) for a detailed discussion of the EOSs including t heir source. The t emperature 



has little effect on the EOS especially in the high pressure regime (jSeager et al.l 120071 ): we ignore 
the temperature dependence of the EOS. This simplifies the equations and their solution, while 
enabling a relatively accurate analysis. 

In this problem, we have five variables: 

1. the iron mass fraction (a); 

2. the silicate mass fraction (/3); 

3. the water-ice mass fraction (7); 
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4. the total mass of the planet (M p ); 

5. and the mean radius of the planet (R p ). 

The variables a, (3, and 7 are not independent of each other. Based on the assumption that the 
planet only consists of iron, silicate, and water we have: a + /3 + 7 = 1, which can also be expressed 
as 7 = 1 — a — p. We therefore have four variables: a,/3,M p and R p . Given any three of these 
variables, we can determine the fourth variable uniquely. We can also see that given an M p and 
R p , there is a relationship between a and (3. There is not a single value of a and (3 that produces 
a given M p and R p . Instead, there are infinite pairs of a and (3 which give the same M p and R p , 
and we call this a degeneracy in the interior composition. 

2.2. Algorithm for Solving the Differential Equations 

Our program integrates from the surface r = R p inward to the center of the planet. The 
outer boundary condition is m(R p ) = M p and P(R P ) = 0. That is, at the surface, the mass is the 
specified total planet mass and the pressure is approximately 0. 

We aim to interpret observations of a planet of a given mass and radius. We therefore choose 
to intergrate inwards instead of outwards based upon the known parameters of the planet (M p , 
R p , and P(R p )). The independent variable is r, decreasing from r = R p to r = 0. The interior 
boundary condition is m{r) = at r = 0. Typically, when integrating m(r), m does not equal zero 
at r = 0. We therefore must iterate, tuning the mass fraction of each layer until m = and r = 
is reached. 

Given a single M p and R p the computer program finds all possible combinations of (a, (3, 
7(= 1 — a — f3)) that give the same planet mass and radius. The program begins with a chosen 
value of a, and then takes a guess of (3. Next, the computer program integrates the differential 
equations (1) and (2) using the above boundary conditions to find the value of the planet radius. 
By comparing this radius to the desired R p , the computer program tunes the value for [3, by using 
the bisection method. This process is repeated several times, until (3 is found to a satisfactory 
accuracy of 1/1000. By varying a within the range of to 1, we can get all possible combinations 
of a and (3 which produce a specific M p and R p . 

2.3. Algorithm for Generating a Database of M p = 0.5-20MQ 

M p and R p are observed parameters and one usually wants to find the corresponding allowed 
a and (3. For a range of M p and R p — corresponding to observational uncertainties — it can be very 
time consuming to use the first algorithm described in £ j2.2i We therefore generate a database that 
is a discrete representation of the relation R p = R p (a, (3, M p ). Figure [Q illustrates the 4D database. 
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This database is a 3-D array which contains the data of R p corresponding to each combination 
of a and (3 (0 to 1 with 1% spacing) and M ffi (ranging from 0.5 to 20 M ffi with O.25M0 spacing). The 
database can be used via linear interpolation to find R p for any given M p , a and (3. A conservative 
estimate of the fractional error in the database interpolation is 1/1000. For the same range of 
M p and R p , interpolation in the database is about 45 times faster than solving the differential 
equations. 

To generate a database of all values of a, (5, and R p for M p ranging from 0.5 to 20 M®. Given 
M p , a, and /3, this algorithm also integrates from the surface inward to r = and m(r) = to find 
R p . In contrast to the first algorithm (which solves for a given M p and R p ), a single integration in 
radius results in the desired solution of R p , for a given M p , a, and f3. In other words is there is no 
iteration required, making this algorithm much more efficient. 



2.4. Instructions for Downloading and Using the Code 

The code is based in MATLAB and can be downloaded from |http : / / web . mit . edu /zengli / www / 1 
under "Research Field" or from Ihttp:/ /seagerexoplanets.mit.edu/research/interiors.html If using 



this computer code please cite this paper and also lSeager et al.l (|2007l ). 



We have made two different codes available. The codes have the same output, but the first is 
based on a differential equation solver ( §2.2[) and the second code is based on interpolation of the 
large database ( §2.3|) . For the codes, the planet mass must be in the range 0.5-20 AT®. The inputs 
to the codes are: the planet mass in Earth masses (M p ), the planet mass uncertainty in Earth 
masses ((Tm p ), the planet radius in Earth radii (R p ), and the planet radius uncertainty in Earth 
radii (<7r p ). The values om p = and (Jr p = are allowed. If the combination of input values M p 
and R p are unphysical, the code will return an error. 

ExoterDE(Af p , <tm p ? Rp, &R P )- This code solves the two differential equations described 
in §2.11 This code consists of three subroutines (each of which must be downloaded) that are 
automatically called by the above command. The first subroutine is the differential equation 
solver, which also reads the equations of state. The second subroutine contains the actual differential 
equations. The t hird subroutin e plots the ternary diagrams; this subroutine calls a ternary diagram 



plotting routine ( Theune 2005 jf| which plots a single line for each of the 1-, 2-, and 3-a contour 



lines. See £ j4.1i An example from this code is shown in Figure HI 

ExoterDB(M p , <tm p ? Rp, &R P )- This code reads in the database of M p , R p and fractional 
composition (a and 0). The output is a ternary diagram, shaded throughout the 1-, 2-, and 3-cx 
contour curves. This subroutine uses the same ternary diagram plotting routine as described above. 



3 http:// www.mathworks.com/matlabcentral/ 
fileexchange/loadFile.do?objectId=7210&objectType=file 
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The differential equation solver ExoterDE is much slower than the database extracter Exo- 
terDB. In principle, ExoterDB is more accurate than ExoterDB. 



3. Data Display 

3.1. 2-D Cartesian Diagram 

For a given M p and R p we want to know the interior composition of the relative mass fraction 
of the three components. There are three variables we have solved for (a, /3, 7), but only two of 
them are independent (since 7 = 1 — a — (3). Therefore points on a 2-D diagram can describe all 
the possible combinations of a, (3, 7 for a given M p and R p . We show such a solution in Figure EJ 
We note that < a < 1, < /? < 1, a + (3 < 1, and therefore not every point in the 2-D plane will 
correspond to a set of a, f3, 7. Only the points which are in a right-angled triangular region will 
respresent the set of allowed solutions. 



3.2. Ternary Diagram 



Ternary diagrams to describe the interior composition of exoplanets were introduced bv lValencia et al. 



(|2007l ). In a ternary diagram, a, (3, and 7 are each one axis of an equilateral triangle. Although 7 
is extraneous, the ternary diagram is useful because it is more intuitive to see the three components 
of the planet interior (in a symmetric way) compared to a 2D Cartesian diagram with only two of 
the components. Figure [3] shows how to read a ternary diagram. 



3.3. Relationship Between the 3D and 2D Cartesian Diagrams and the Ternary 

Diagram 

To explain the full origin of a curve on the ternary diagram we start with the 3D Cartesian 
diagram with all solutions of M p , R p and composition (in terms of a and /3), as shown in Figure [TJ 
We take an isoradius and isomass surface as shown in Figured! As an example, in Figured] the red 
surface is the iso-radius surface of R p = 1.7 and one of the blue colored planess is the iso-mass 
surface of M p = 4M ffi . These two surfaces intersect each other and result in a curve. This curve can 
be projected vertically to the x-y plane which is the iso-mass plane. We therefore have a (iso-radius 
and isomass) curve on the isomass plane. This curve is shown in Figure [2] in a Cartesian diagram. 

The Cartesian and ternary diagrams are two different ways to represent the same information. 
There exists a linear coordinate transformation between the two. That means if a function is a 
straight line appearing in the 2-D Cartesian diagram, it will still be a straight line in the ternary 
diagram. 
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The transformation from 2-D Cartesian coordinates to the ternary diagram coordinates is 

^ternary — 2 ^ 
V ternary — <y (1 X ?/). (4) 

Here x and y are the coordinates of a point in a 2D Cartesian diagram, the xt ernary and ytemary 
are the coordinates of the point in ternary diagram in the Cartesian grid variables. Figure [3] shows 
the same R p = 1.7 Rq, M v = 4M® curve represented by a ternary diagram. 



4. Results and Discussion 
4.1. Observational Uncertainties 

Real planet mass and radius measurements have uncertainti es. The planet mass and radius 



uncertainties are typically 5 to 10 percent (e.g., ISelsis et al.l 120071 ) . and even smaller for the most 



favorable targets. We now present examples of ternary diagrams that include the mass and radius 
uncertainties. 

We consider uncertainties of 1, 2, and 3 standard deviations (a) from the measured value. See 
Figures d] and [5j In more detail, the uncertainty in composition on the ternary diagram is 



\J (PcompM) 2 + {VcompR) 2 , (5) 



where comp refers to composition and compM and compR refer to composition uncertainties caused 
by the planet mass and radius uncertainty respectively. Here we have assumed that the uncertain- 
ties in mass and radius are independent from each other and have assumed the linearity of the 
superposition of small uncertainties. 

Figure 0] shows a planet with M p = 10 ± 0.5M© and R p = 2 ± 0.1Af e . We see that taking 
the 3-<T limit, almost the entire ternary diagram is filled. In other words, for a 5% 3-er (i.e., 
15%) uncertainty on the planet mass and radius, the interior composition in terms of fractional 
composition of iron, silicates and water cannot be determined. The reason this example fills the 
almost the whole ternary diagram is that a planet with 10M® and 2R^ has an average density in 
between two extreme cases (purely iron or purely water). Therefore, a large variety of different 
combinations of iron, silicate and water can result in a similar M p and R p . Even taking a 1-a 
uncertainty of the planet mass and radius, the uncertainty in internal composition is large. 

We note that an uncertainty in R p has more of an effect on the uncertainty in the interior 
composition than an uncertainty in M p . This is because the planet's average density p ~ M p /R s . 
Considering error propagation, the uncertainty in radius has a three times larger effect on the 
uncertainty in average density than does the mass uncertainty. 
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We show ternary diagrams for planets with various masses, radii, and 5% fractional uncertainty 
in Figure [5j Only solutions in part of the ternary diagram are allowed, despite considering the 3-<r 
range. In Figure [5]b, the upper 3-cr boundary is absent because it goes below the lowest allowed 
density of the 2M ffi and 1.5-R® planet and is thus unphysical. In Figure [5fc we see the opposite case, 
where the lower 3-cr boundary is absent because it goes above the highest planet density allowed 
and is thus unphysical. 

A ternary diagram for a fixed planet mass and radius that includes observational uncertainties 
is one of the primary outcomes of this paper. 



4.2. Model Uncertainties 

The model and computer code we present assumes a differentiated planet composed of an iron 
core, a silicate mantle, and a water ice outer layer. The division into three major materials is 
based on the point that the densities of iron, silicate, and water are much more different from each 
other than any minor compositional variant of e ach individua l mate rial. The model neglects phase 



variation and temperatures which, as argued in lSeager et al.1 ( 20071 ). have little effect on the total 



planet radius (to an uncertainty of about ~1 to 3 percent uncertainty in planet radius, decreasing 
with increasing planet mass.). 

Low-pressure phase changes (at < 10 GPa) are not important for a planet's radius because for 
plausible planet compositions most of the mass is at high pressure. For high pressure phase changes 
we expect the associated correction to the equation of state (and hence derived planet radii) to 
be small because at high pressure the importance of chemical bonding patterns to the equation of 
state drops. 

Regarding temperature, at low pressures (< 10 GPa) in the outer planetary layers, the crystal 
lattice structure dominates the material's density and the thermal vibration contribution to the 
density are small in comparison. At high pressures the thermal pressure contribution to the EOS 
is small because the close-packed nature of the materials prevents structural changes from thermal 
pressure contributions. 

Although the code can model radii for planets in the mass range .5 to 20 M®, the model is 



more accurate for planets above a few Earth masses I S eager et al.l (|2007l ) 



The model also neglects variation in composition, such as a light element in the iron core as 
Earth and Mercury are believed to have. The model also omits other impurities in the mantle and 
water layer, including iron in the mantle. Molten cores have also been omitted. At the present 
time, these model uncertainties are expected to have an effect on the planet radius much less than 
the ~5 percent radius observational uncertainty. 



For all of the above reasons, we therefore argue that for the present time the observational un- 
certainties dominate the model uncertainties; the model presented here is adequate for an estimate 
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of planet bulk composition. In any event, the main results of our work described in the following 
subsections despite any model uncertainties. 

It is possible to rule out parts of the ternary diagram as being physically unplausible (e.g., 



Valencia et al.ll2007l ). This is based on the initial composition of the protoplanetary nebula and on 
planet differentiation. For example, a pure iron planet is unlikely to exist, because removing all 
of the mantle would be difficult. A pure water planet is also unlikely to exist. Where water ice 
forms, so do silicate-rich and iron-rich materials, making planet accretion of pure water unlikely. 
We prefer to leave the omission of parts of the ternary diagram to users of the code, because in 
exoplanets surprising exceptions to the "rules" of planet characteristics are not uncommon. 



4.3. Spacing, Shape, Direction, and Rotation of Curves on the Ternary Diagram. 

We now turn to a discussion of the spacing, direction, and shape of the curves in the ternary 
diagrams. A quantitative and qualitative description of these is a main point of this paper. We 
emphasize that each curve shown in our Figure [5] diagrams represents a different mass and radius. 
The curves to the lower right are more dense, as they have a higher mass and lower radius than 
the curves moving to the upper left. 

All behavior results from the equations of state of the materials, and, in some cases, how they 
behave differently under pressure. 

We begin with an equation that we use repeatedly in this section. We consider the simplified 
case that the planet core has a uniform density, where pp e is the average density of Fe in the core, 
PMgSi0 3 is the average density of silicate in the mantle, and pn 2 o is the average density of water 
in the outer water layer. We then have 



3^ 

where a. + /3 + 7 = 1. 



aM p | (3M p | 7 M P 

PFe PMgSiOz PH 2 



(6) 



4-3.1. Spacing 

The iso-mass iso-radius curves for adjacent curves with equal differences of mass and radius 
have uneven spacing on a ternary diagram. This spacing is generally smaller in the lower right 
part of the ternary diagram (high iron fraction region) than in the upper left part (high silicate 
or water fraction region). The curves in the lower right part of the diagram have a higher density 
(higher mass and smaller radius) than the curves found on the upper left part of the diagram. The 
density is not a linear function of both mass and radius hence we do not expect equal spacing on 
the ternary diagram. We can, however, give both a quantitative and qualitative explanation of the 
uneven spacing. 
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We can provide a quantitative description, beginning with equation ©, but using a = 1 — (3 — 7 
to get 



(1-/9-7) 

PFe 



ft 



1 



+ + 

PMgSi0 3 PH 2 



(7) 



We use the fact that the curves on the ternary diagram are almost perpendicular to the water 
side of the ternary diagram, and therefore set the silicate mass fraction ft = for our discussion. 
In other words, the distance (separation) between the points produced by the intersection of the 
isomass-isoradius curves and the water axis is a good representation of the spacing between the 
curves throughout the ternary diagram, 



1 f 1 1 

— + 7 _ — 

PFe \PH2O PFe 



(8) 



We define the following constants: 



At 



PFe 



and 



Equation ([8]) then becomes 



We note that A\ > and A 2 > 0. 



An 



PH2O PFe 



4/3irRl 



(9) 
(10) 

(11) 



Now we proceed to take the derivative of equation (jlip with respect to the change of mass 
(dM p ) and change of radius (dR p ). 



4"7T 

^2^7 = Y 



3 1^ tl 1^ 

j£ dR >-iil dM > 



(12) 



We can also rewrite equation (jlip in terms of the overall average density of the planet (p) , 

7A2 = --Ai, 
P 

and the corresponding derivative relative to the overall average density (p) 



(13) 



A 2 d^y = - —zdp. 
P 2 



(14) 



This leads to our quantitative understanding, where we first recall that dj is the water fraction 
spacing on the ternary diagram. The lower right part (Fe-rich) of the ternary diagram is where 
the average density (p) of a planet is high. This implies (for the same dp) dj is small (since p is in 
the denominator). In the upper left region of the ternary diagram (water-rich) the average planet 
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density is smaller than a planet located in the lower right part of the diagram, and therefore d'y is 
larger. 

Qualitatively, to have wider spacing an increasing water fraction is needed. In other words, 
towards the upper right part of the ternary diagram, for the same density difference more water 
than iron must be replaced. 

We note that the spacing is predominantly the result of the nonlinearity of equations (jlip and 
(I14h and the average density of each compositional layer, not of any P-p properties of the EOS 
(i.e., how materials condense under high pressure). This statement is correct under the assumption 
of a single material for each layer in the planet (in our case for iron, silicate, and water). The 
assumption that the average density within each layer does not change significantly from curve to 
curve (for example, from the 1-a curve to the 3-a curve for the case M p = 2M^ and R p = 1.5Rq) 
is reasonable for a ternary diagram that spans only a small mass and radius range. 



4-3.2. Shape, Direction, and Rotation 



To explain the shape and direction of the curves in the ternary diagrams we start by explaining 
the slope of the curves in the Cartesian diagram. In other words, we are aiming for an expression 
of da/ dp. 



We start with a different form of equation 



a 



1 



1 



PH 2 PFe 



+ P 



1 



1 



PH 2 PMgSiQ 3 



1 



AirR* 
Ph 2 o 3M p 

We now differentiate this equation with respect to a and (3 to find 



da 



1 



1 



PH 2 PFe 



+ df3 



1 



1 



PH 2 Q PMgSiQ 3 



0, 



(15) 



(16) 



based on the assumption that the average density is changing slowly with respect to the change in 
composition. 

We write the slope of the curves in the Cartesian diagram 



d/3 



da 



Ph 2 o 

PFe 



Ph 2 o 



We can see that 



dp 
da 



PMgSi0 3 



> 1, 



(17) 



(18) 
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assuming that p H2 o < PM g SiO a < PFe 



It can be shown that the slope in the the ternary diagram has a positive connection to the 
slope in the Cartesian diagram. Using the the equations that convert the Cartesian coordinates to 
coordinates on the ternary diagram (equations ([3j) and (Jl|)) we find that 



dp, 



ternary 



da t 



ernary 



+ 1 



(19) 



We can now go on to describe the direction, shape, and rotation in the Cartesian diagram 
based on equation (|17p . with the knowledge that the same qualitative behavior will appear in the 
ternary diagrams. We first emphasize that the slope of an isomass, isoradius curve on the ternary 
diagram describes adding and removing mass of the different species. 

We begin with a qualitative explanation of the direction of the curves on the ternary diagram — 
why each curve goes from the lower left to the upper right. This is largely a coincidence in the 
different values of ph 2 o, P~MgSi0 3 , and pFe- The coincidence lies in the fact that iron is more dense 
and water less dense than silicate — and for zero-pressure densities, an equal mass of iron and water 
combined densities are roughly similar to the silicate density. At the lower left of the ternary 
diagram in Figure O the silicate mass fraction is 80 percent, the iron mass fraction is zero and 
the water mass fraction is 20 percent. As the silicate fraction decreases, a combination of equal 
parts iron and water must be added to maintain the same overall planet mass and radius. This 
description is consistent with the direction of the curves in the ternary diagrams. 

More quantitatively, from equation (fTT|) . we can take the zero-pressure densities of Fe, H2O, 
and MgSi03 to find the slope of the curve on the cartesian diagram 



and from equations ([3|) and (jlj), 



d£_ 

da 



dPternary 



1 



(20) 



(21) 



dOit ernar y 

We can consider removing a fixed amount of iron mass, e.g., 1 gram. According to equation (|20p . 
2 grams of silicate must be added. For mass balance (because each curve on the ternary diagram 
represents a planet of fixed mass and fixed radius), 1 gram of water must be removed. The direction 
of the curves on the ternary diagram do correspond to removing roughly equal masses of iron and 
water for every mass of silicate added. 

Rotation 

In Figure (5^,-d (where a through d are in order of increasing mass), we see that for more 
massive planets, the isomass-isoradius curves are rotated. In other words, the slope of the curve 
on the ternary diagram increases for increasing mass. 
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We again return to equation (fT7|) . With increasing planet mass, each of Ph 2 Oi PMgSiOz-, and 
pp e changes. Because iron is in the core, it suffers more compression than the silicate mantle or 
water icy layer. In other words, pFe must increase more than pH 2 o and pMgSiO s as the planet mass 
increases. Therefore, the numerator of equation (|17p gets larger and thus the slope of the curve 
increases. The rotation is counter clockwise with increase of mass. 

Slope of the Curves 

We now turn to discuss the slope of each isomass-isoradius curve. We see from Figures [4] and 
[5] that the slope of each isomass-isoradius curve for a fixed mass and radius is greater in the lower 
left part of the curve than in the upper right part of the curve. 

The slope is due to the differential compression of water, silicate and iron under pressure. 
This slope is again explained by equation (|17p . At the lower left the slope is smaller; this is the 
silicate-rich region of the ternary diagram. There is more silicate and less iron and water. The 
PMgSiOs will increase as it gets compressed in the inner part of the planet. This causes the slope 
in equation (|17h to get smaller. In contrast, at the upper right, there is little silicate, but more 
water and more iron. The silicate and water are less compressed, but the iron is more compressed, 
making the slope increase. 

For a conceptual explanation, first recall the idea described above that removing 2 grams of 
silicate can be compensated by adding approximately 1 gram of water and 1 gram of iron. In a 
silicate-rich planet (lower left of the ternary diagram), silicate is compressed. For a massive planet, 
this compression makes the silicate density closer to iron's density than in the uncompressed case. 
Therefore, removing a fixed mass of silicate requires much more iron than water to be added, for a 
fixed mass and radius. 

In contrast, along the upper right part of the ternary diagram, the planet is iron-rich or water- 
rich. The slope of an isomass-isoradius curve is steeper than a curve in the lower left part of 
the ternary diagram. For an isomass-isoradius curve in the upper part of the ternary diagram, the 
planet has more iron, the iron is very compressed (There is less silicate and the silicate is overall less 
compressed compared to a planet in the lower left of the ternary diagram) . If a fixed mass of silicate 
is removed, more water than iron must be added to compensate for the density of compressed iron. 
The density of water does not change much, because water is always the outer layer and thus the 
least compressed. 

As a qualitative explanation, if the average density of each of the three layers remains constant, 
then an isomass-isoradius curve should always be a straight line in either the Cartesian or ternary 
representation. The curvature in a given isomass-isoradius curve appears because of the compression 
of material under pressure, which is a property of the EOS. 
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4.4. Shape and Direction of an Isoradius Surface in our 3D Representation 

We return to the 3D representation of the relationship between mass, radius and composition 
shown in Figure [TJ An isoradius surface (red oblique surface) is shown for R p = 1.7 R®. An isomass 
surface (one of the blue flat planes) is shown for M p = -iM^. At the bottom tip of the isoradius 
surface, the silicate mass fraction = iron mass fraction = 0, and the planet is composed of 100% 
water. This is the minimum mass for this radius. The isoradius surface also has a maximum mass, 
reached by a composition of 100 percent iron. 

For the same radius, if either the silicate or iron mass fraction is increased, the planet mass 
must also increase. This is because iron and silicate are denser than water ice. We further note 
that for an increase in iron mass fraction, the mass of the planet must increase more steeply than 
for an increase in the silicate mass fraction. This is seen by the different length and shape "edges" 
of the isomass surface in the iron-mass plane and silicate-mass plane in Figure [TJ 

We now show that the shape of the isoradius surface is concave. The isoradius surface can be 
considered as the isovolume surface, where the volume is the sum of the core, the silicate mantle, and 
the water crust. To calculate the total mass of a point on the isoradius curve, we use equation ©, 
and since R p is constant on this surface we can rewrite this equation (with C a constant) as, 



M p = C 



<* | P | (l-a-P) 



_PFe PMgSiO-s PH 2 

The term in brackets is a linear function and its inverse is hyperbolic. For example, if we let a 
(no iron), then we have 

-l 



(22) 



M p = C 



P 



1 1 



(23) 



_PH 2 \PH 2 PMgSiOz. 

where 

C = -ttR^. (24) 
Both terms in the square brackets are positive and < (3 < 1. This is in the form 

where 

Ci = ^— (26) 

PH 2 

and 

C 2 = -_ (27) 

PH 2 PMgSi0 3 

We have C\ > C 2 > 0. This is the form of a concave hyperbola, because, taking the first derivative 
we find 

dM p _ CC 2 (28) 



dP (Ci - C 2 P) 



2 
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The slope increases as (3 increases because C\ — decreases but is always greater than zero. 

Although we have made the simplification that the average densities of each layer remain 
constant, our qualitative description holds because pn 2 o < PMgSi0 3 < PFe as long as the planet is 
differentiated into layers of increasing density towards the planet center. 

A similar argument shows that the total iron fraction vs. total mass is also a hyperbola, 
making the whole isoradius surface concave. 



5. Summary and Conclusion 

An ambiguity in an exoplanet interior composition remains for any planet with a measured 
mass and radius, no matter how precisely measured. W e can accept this ambiguity and quantify it 



with the aid of ternary diagrams (j Valencia et al.1 120071 ). We have presented ternary diagrams for a 



single planet of fixed mass and radius, for a planet composed of an iron core, a silicate mantle, and 
a water ice outer layer. Our ternary diagram presentation includes observational uncertainties. We 
have provided a publically available computer code to generate a ternary diagram for a given input 
mass, radius, and observational uncertainties. 

In addition to presenting ternary diagrams for fixed mass and radius, we showed their origin 
from a 4D database (M p , R p , iron mass-fraction a, and the silicate mass fraction (3; recall that 
the water-ice mass fraction 7 = 1 — a — /3). We further described the shape and direction of the 
composition curves on a ternary diagram . 

We conclude with the sentiment that in order to fully understand the interior structure of an 
exoplanet, a third measurement beyond planet mass and radius is required. 
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Fig. 1. — The 3D representation of planet mass, planet radius, and composition expressed as mass 
fractions of iron (a) and silicate (/?). The x-axis is the iron mass fraction and the y-axis is the 
silicate mass fraction. The water mass fraction comes from a + (3 + 7 = 1. The z-axis is the 
total planet mass in Earth masses. The planet radius is not indicated, except by three separate 
surfaces of constant radius. An isomass and isoradius surface intersect with a curve in a 2D iron 
mass- fraction vs. silicate mass fraction Cartesian diagram (Figure [2]). 
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Fig. 2. — The Cartesian diagram for silicate mass fraction (/?) vs. iron mass fraction (a) for a planet 
with M p = 4M ffi , R p = 1.7-R© (Note that the water mass fraction 7 = 0.4 since 7 = 1 — a — (3.) 
A planet of a fixed internal composition is represented as a point in the grey shaded region (only 
compositions in the grey shaded region are allowed.) With a given planet mass and radius there 
is not a unique interior composition; the allowed compositions are described by the curve. The 
curve originates from the intersection of the isomass surface and the isoradius surface in the 3D 
representation (Figure [1]). 
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Fig. 3. — A ternary diagram for a planet of the same mass and radius as shown in the corresponding 
Cartesian diagram in Figure [2j The light grey x and y axes are shown to illustrate how this ternary 
diagram relates to a Cartesian diagram; xt ernar y and ytemary are the coordinates transformed from 
the Cartesian coordinates. (See § 13,31 for the conversion equations.) Any point on the curve is a 
possible combination of iron, silicate, and water that will result in a planet with with M p = 4M ffi 
and R p = 1.7R®. This Figure also illustrates how to read a ternary diagram. Consider a triangle 
oriented such that the value 1 (of a given material) is at the triangle's apex. The mass fraction of 
that given material can be read off of a horizontal line that is perpendicular to the line connecting 
the triangle's apex with the triangle base. 
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a (iron) 

Fig. 4. — A ternary diagram including the mass and radius uncertainties for a planet of a fixed 
mass and radius. This example is for a planet with M p = 10±0.5M ffi and R p = 2±0.1-M®, showing 
the 1-, 2- and 3-a uncertainty curves. This Figure shows the continuous distribution of possible 
combinations of iron, silicate, and water with increasing uncertainties according to the color bar. 
Notice that considering the 3-cr uncertainties almost the entire ternary diagram is covered — in other 
words there is no constraint on the planet internal composition. See text for a discussion of the 
direction and spacing of the curves. 



- 20 - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



a (iron) a (iron) 

Fig. 5. — Ternary diagrams including 5% mass and radius uncertainties for planets of fixed mass 
and radius. Panel a) M p = 1 ± 0.05M® and R p = 1 ± 0.05M®. Panel b) M p = 2 ± 0.1M® and 
R p = 1.5±0.075M©. Panel c) M p = 8±0.4M ffi and R p = 1.5±0.075Af©. Panel d) M p = 16±0.8M e 
and R p = 2.5 ± 0.125M ffi . See text for a discussion of the direction and spacing of the curves. 



